Instabilities of wave function monopoles in Bose-Einstein condensates 



Th. Busch, J. R. Anglin 
Institut fiir Theoretische Physik, Universitdt Innsbruck, A-6020 Innsbruck, AUSTRIA 



00 

> 

o 



> 



00 
On 



o 
o 



We present analytic and numerical results for a class of 
monopole solutions to the two- component Gross-Pitaevski 
equation for a two- species Bose condensate in an effectively 
two-dimensional trap. We exhibit dynamical instabilities in- 
volving vortex production as one species pours through an- 
other, from which we conclude that the sub-optical sharpness 
of potentials exerted by matter waves makes condensates ideal 
tools for manipulating condensates. We also show that there 
are two equally valid but drastically different hydrodynamic 
descriptions of a two-component condensate, and illustrate 
how different phenomena may appear simpler in each. 

PACS number (s): 03.75.Fi, 03.65 Ge 



Monopoles are particle-like field configurations with 
which one can associate a topological charge. As per- 
haps the most obvious way of making a collective par- 
ticle out of a condensate field [|l|-^, they are a natural 
goal for condensate state engineering. In this Letter we 
discuss a specific case D = 2 of a class of /^-dimensional 
monopoles, previously introduced in general and exam- 
ined for D = 1 Although the notion of topological 
charge (in D = 2, the 'winding number') is useful in iden- 
tifying and classifying these particle-like configurations, 
in our case this charge is not guaranteed to be conserved. 
Indeed the structures we discuss are quite prone to dy- 
namical instabilities. They are nevertheless also robust, 
in that their instabilities typically do not destroy them, 
but rather distort them, in ways which are both intrin- 
sically interesting, and illustrative of some basic features 
of multi-component superfluid dynamics. 

Our paper is organized as follows. After presenting 
our concept of a 'wave function monopole', we consider a 
one-dimensional 'ring monopole' which bears much the 
same useful analogy to the 2D monopole as the one- 
dimensional persistent current does to the vortex (and 
which also merits serious consideration in its own right). 
In this context we briefly discuss the consequences of hav- 
ing different scattering lengths among the two atomic 
species. We then present and analyze the results of nu- 
merical solutions to the two-component Gross-Pitaevski 
equation. The specific cases we examine are chosen, 
mainly for theoretical simplicity, from among a vast 
range of possibilities, and while we have also investigated 
enough other cases to confirm that the results we present 
are not qualitatively altered by slight changes in param- 
eters, we are not insisting that our precise examples will 
be optimal laboratory subjects in the near future. We do 
identify some basic features of a class of structures that 
are realistic targets for experiments, and so our paper 
concludes with a few remarks of experimental relevance. 



The basic concept of a monopole in D dimensions is 
sketched in Fig. 1 (for D = 2). The arrows represent an 
order parameter having as many 
I ^ components as the dimensional- 

ity of the space in which the 
^ monopole exists. There is much 

~ freedom in the actual identifica- 

// y tion of this order parameter. For 

^ /i ' ^ ' instance, in two dimensions its 

• ' magnitude and direction could 

FIG. 1. Monopole. represent the modulus and phase 
of the macroscopic wave function 
of a one-species condensate; this 'monopole' is the well- 
known superfluid vortex. In the case we now consider, 
however, the two dimensions of the arrows instead rep- 
resent the wave functions of each component of a two- 
species condensate, in a configuration in which both com- 
ponents happen to be real: ('^i,'02) = /(r)(cos ^, sin^), 
where r, are the usual polar co-ordinates . 

As with the vortex, the modulus of the order parameter 
is the density of the condensate (here, the total density 
of both components), so that the 'particle-like' core of 
the monopole is in fact a local minimum of density: a 
void or bubble, maintained by destructive interference of 
matter waves. As with the vortex, however, it is instruc- 
tive first to avoid the core r ^ of the monopole, and 
consider only its behaviour at large r, where the radial 
dimension becomes unimportant, leaving the effectively 
one-dimensional problem of a two-species condensate on 
a circle of fixed radius R. So we replace /(r) 1, and 
examine the 'ring monopole' (V^i,'02) as a stationary so- 
lution to the one-dimensional Gross-Pitaevski equation 
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where /i = 1 — (2gR^)~^ provides = 0. 

For the ring monopole we can obtain analytically 
the Bogoliubov spectrum of perturbations, ijjj ipj + 
e(j){0^t). Working to linear order in e yields the modes 
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The ± index distinguishes acoustic and optical branches. 
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in which the two species' density perturbations are re- 
spectively in and out of phase: 
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Note that rotating these modes by ^ ^ (^+7r/2) produces 
an independent set of modes, with the same frequencies; 
every frequency is thus two- fold degenerate. 

The modes k = (for which the ± branches co-incide) 
and k± = 2— provide the four zero modes due to the 
model's U(2) symmetry (which already includes spatial 
rotation). All the other modes are positive frequency 
excitations, except for the two k± = 1— modes; and 
for gR^ > 3/4, these become dynamical instabilities. In- 
finitesimal excitation of these modes then grows exponen- 
tially until it becomes finite, and the Bogoliubov theory 
is inadequate. To follow the evolution into this regime, 
we solve the Gross-Pitaevski equation numerically, using 
the split operator method as described in [^. It turns 
out that the finite perturbation does not mix with other 
modes to produce irreversible (long revival time) decay, 
but grows to a maximum size, then shrinks back; and 
the cycle repeats. This indicates that the two unstable 
modes form an isolated subsystem of two degrees of free- 
dom, with a 'Mexican hat' potential. 
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FIG. 2. Dynamical instability of the ring monopole, for 
gR^ = SOO/tt. Five times are shown, the configurations at 
the last two times being precise revivals of those at the first 
two. The top three plots show |'0i|^ (broken) and |'02|^ (solid); 
the middle three show IV^+I^ and for the same configu- 

rations. (In the left middle plot the small initial perturbation 
can just be seen.) The three rings at bottom are overhead 
views of the iV^i |^ + |'02|^, with darker shading for lower total 
density. The location of the dark spot is determined by the 
initial conditions, but if a small random potential is added on 
the ring, over many cycles the system performs what might 
be called 'superfiuid roulette'. 



Before going on to examine instabilities of monopoles 
in two dimensions, we pause here to address an issue 
of generalization, which will lead to the development of 
some concepts useful for understanding two-dimensional 
phenomena. In (^ it is assumed that the three scattering 
lengths and two particle numbers in the problem are all 
equal; this is very close to the case for ^^Rb In fact 
the simple solution we examine here can be generalized to 
encompass all positive scatteringlengths gij. When sep- 
aration is favoured (det ^ < 1) one can obtain solu- 
tions in which our cosine and sine are replaced by Jacobi 
elliptic functions. (Phase-separating monopoles in two 
dimensions resemble crossed domain walls, dividing two 
condensates each with a Josephson 7r-j unction through 
the monopole core. They appear to have qualitatively 
similar behaviour to the simplest case with gij = g.) 
When mixing is favoured (det ^ > l)^ o^^^ can also gen- 
eralize monopoles, but in a somewhat surprising way. If 
one merely changes basis, defining = -^{i^i i ^'^2), 
one discovers that our monopole can also be described 
as two superimposed, counter-rotating vortices (SCV). 
Ghanging scattering lengths so as to favour mixing of 
the two species will then stabilize this configuration. 

Of course, this change of basis is only a symmetry of 
the system when all scattering lengths are equal (and 
hence only then is it really true that the monopole and 
the SGV are the same object). Nevertheless the two 
bases can always be used, and so one has two very dif- 
ferent hydrodynamic pictures which describe the same 
superfiuid physics. This is our first principle, that with 
two-component superfiuids one has two equally valid hy- 
drodynamic descriptions, and that one is free to base 
intuition on either of them. As will be apparent below, 
the differences between them can be dramatic. 

We now return to the case where all scattering lengths 
are equal, and investigate the monopole in two dimen- 
sions. We note that ('0i,'02) = /(r)(cos^, sin^) is a sta- 
tionary solution to the Gross-Pitaevski equation (GPE) 



as long as /(r) satisfies the same nonlinear equation as 
the modulus of the two-dimensional vortex in the trap: 
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)/ = 2g{p - y)f . 



(5) 



We have studied such monopoles in harmonic traps 
numerically, using the split-operator method to solve the 
two-component GPE in two dimensions, in both real and 
imaginary time. We have examined a wide range of differ- 
ent parameter regimes, and compared results with grids 
of varying sizes; and while below we show displacements 
and motions which are exactly in one of our grid direc- 
tions, we have also checked that the same behaviour oc- 
curs at arbitrary angles. 

If a monopole is displaced from the centre of the trap, 
it begins to fall back towards the centre; but this does 
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not mean that it is stable. By considering the hnearized 
theory in imaginary time, one can show that there are 
lower energy configurations near the monopole, reached 
by short range perturbations near the core. While this 
analysis does not suffice to determine whether the insta- 
bility is dynamical or merely energetic, numerical study 
of small amplitude oscillations of a monopole in a trap 
indicates that monopoles are dynamically unstable: the 
amplitude of small oscillations grows exponentially. 
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FIG. 3. Total density p = iV^il^ + |'02|^ in the centre of a 
harmonic trap, for a monopole initially displaced very slightly 
from the centre, at g J(Px\i/jj\'^ = 2000. Note logarithmic 
vertical scale. Successive minima do not reach p = 0, because 
growing speed of the monopole involves filling in of the core. 

The instability, when infinitesimal, is a self- amplifying 
oscillation of the monopole back and forth through the 
centre of the trap; this is at least qualitatively similar to 
the instability on the ring. As the motion grows finite, 
vortex pairs form. As can be seen from Fig. 4, the motion 
of the monopole involves the pouring of one of the two 
species through a narrow channel formed by the other 
species (which is comparatively passive throughout). As 
shown in Fig. 5, very rapid flow develops along the sides 
of this channel, until vortex pairs suddenly nucleate at 
the maxima of the passive species' density. Each pair 
then separates, into an inner vortex pinned by the pas- 
sive species, and an outer vortex, which gradually drifts 
away. Fig. 6 shows this process as the growth of a pair 
of branch cuts in the condensate phase. It may be quite 
hard to directly detect vortex pairs formed in this man- 
ner; but their appearance vividly illustrates a potentially 
useful fact: that sharp-edged potentials (such as most 
readily generate vortices) are much more easily formed 
with matter waves than with light waves. If one wishes 
to manipulate condensates with high resolution forces, 
therefore, the best tools are probably other condensates. 

Despite spawning vortices, an initially displaced 
monopole continues its fall towards the centre of the trap 
with little apparent impediment. As (or just before) it 
reaches the centre, however, its motion changes. In most 
cases we have investigated, it abruptly slows (though in 
some cases, and for reasons which are not clear to us, it 
seems instead to accelerate, and then decay rather vio- 
lently). Thereafter, the minimum in total density at the 
monopole core gradually fills in, until only a slight depres- 
sion remains. But even at this late stage, a double-lobed 
pattern of the two densities (as in Fig. 4) remains; and 



in each species a phase difference of tt persists across the 
monopole, although the sharp jump in phase has been 
smoothed out. The instability appears to have satu- 
rated, leaving a 'smeared' monopole which appears to 
be dynamically (although not energetically) stable. This 
eventual stability is much easier to understand in the 
SCV basis. There, all that happens is that the super- 
imposed counter-rotating vortices slip apart, becoming 
a pair of reciprocally pinned vortices: each condensate 
species has one vortex, and each also fills a strong po- 
tential well formed by the other species' vortex. (This 
may be seen as an analogue of the SCV instability on the 
ring.) Once this separation has occurred, the stability of 
two pinned vortices is unsurprising. See Fig. 7. 
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FIG. 4. iV^il^ (top) and |'02|^ for a monopole initially 
displaced in the direction ^ = tt, at early and late times 
(g Jd'^x\tl;j\'^ = 2000). Note 'throttle' formed by |^2r 
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FIG. 5. Superfluid velocity component Vxi = ^(^^InV^i), 
at various times in the same simulation illustrated in Fig. 4. 
Plots are 'edge-on' views of the velocity field Vxi{x, y), looking 
in the positive x direction: the velocity is very small every- 
where except near the arc on which vanishes, so the 
edge-on view is convenient.) Initial fluid velocity is negative 
as the monopole moves in the positive {0 = 0) direction. 
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FIG. 6. Phase of ipi{x,y) before and after appearance of 
vortex pairs. The branch cuts are indeed jumps of 27r. 
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FIG. 7. Monopole instability in the SCV basis. Top fig- 
ures show total density p — IV^+I^ + IV^-I^ at early and late 
times; bottom figures show IV^+I^ (IV^-I^ being the mirror im- 
age in the x-axis.) No new vortices form; instead, the initially 
superimposed vortices slip apart, and each is pinned. 

On the other hand, the monopole basis can be more 
illuminating of other phenomena. Two directly aligned 
monopoles perform a curious dance: they attract each 
other, merging into a single monopole of winding number 
two; this then separates again, but in the perpendicular 
direction; and the cycle repeats. In the SCV basis, this 
is a rather confusing system of orbiting vortices, and pin- 
ning peaks that form and dissolve; but in the monopole 
basis, it is quite simple (see Fig. 8). And the peculiar 
scattering of monopoles at right angles in a head-on col- 
lision is actually what is expected for monopoles in gauge 
theories as well We find analogous dances, with 7r/n 
scattering, for monopoles of higher winding number n. 



the two species can be imaged separately, the pattern of 
density lobes is large and obvious, and clearly indicates 
the presence of a monopole core. And finally, the two 
principles of alternative hydrodynamic pictures, and us- 
ing condensates to manipulate condensates, may be of 
practical value in a wide range of future experiments on 
multi-component condensates. 
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FIG. 9. Adiabatic transfer of population (see Ref. [10]), 
from an internal atomic ground state and the motional ground 
state in the trap, to the monopole configuration of two excited 
internal states. We take Aq/cl = 0.1, g J d^x\ipjf — 1000, de- 
tuning swept linearly from -1 to +5 trap frequencies. Further 
parameter optimization may allow total transfer, but there 
may also be advantages, for stabilization and observation, to 
leaving a small ground state population in the core. 
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FIG. 8. Two aligned monopoles scattering in a harmonic 
trap. Top figures show total density; bottom figures show 
iV^il^ and |'02|^ separately. Since |'02|^ essentially does not 
change, it is merely indicated as uniformly shaded lobes. 

We conclude our overture to monopoles in Bose con- 
densates with three remarks on experimental possibili- 
ties. First, despite their instability, creation of structures 
like these, in spherical or toroidal traps, may neverthe- 
less be possible using the adiabatic passage technique of 
Dum et al. See Fig. 9. Second, unlike vortices, de- 
tection of monopoles should not pose a problem: since 
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